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Our modified gravity theory (MOG) was used successfully in the past to explain a range of 
astronomical and cosmological observations, including galaxy rotation curves, the CMB acoustic 
peaks, and the galaxy mass power spectrum. MOG was also used successfully to explain the unusual 
features of the Bullet Cluster 1E0657-558 without exotic dark matter. In the present work, we derive 
the relativistic equations of motion in the spherically symmetric field of a point source in MOG and, 
in particular, we derive equations for light bending and lensing. Our results also have broader 
applications in the case of extended distributions of matter, and they can be used to validate the 
Bullet Cluster results and provide a possible explanation for the merging clusters in Abell 520. 

PACS numbers: 04.20.Cv,04.50.Kd,04.80.Cc,45.20.D-,45.50.-j,98.80.-k 



I. INTRODUCTION 

Our modified gravity theory (MOG), also known as 
Scalar-Tensor- Vector Gravity (STVG, [l|), is a theory 
based on the action principle. The action incorporates, 
in addition to the usual Einstein-Hilbert term associated 
with the metric 5^^, a massive vector field 0^, and three 
scalar fields representing running values of the gravita- 
tional constant G, the vector field's mass and its cou- 
pling strength uj. The vector field is associated with a 
fifth force charge that is proportional to mass-energy. 
This fifth force is repulsive; for a large source mass, at 
large distances, gravity is stronger than that predicted 
by Newton or Einstein, but at short range, this stronger 
gravitational attraction is canceled by the fifth force field, 
leaving only Newtonian gravity. The theory has been 
used successfully to account for the rotation curves of 
galaxies 0, H, 0, [1| , the mass profiles of galaxy clusters 
, and cosmological observations 0, Q without exotic 
dark matter. 

The theory was also used to offer an explanation for 
the unusual features of the Bullet Cluster 1E0657-558 
Q.In 0, a lensing formula was used that was deduced 
from the nonrelativistic acceleration law for test particles 
in the vicinity of a MOG point source. 

In the present work, we derive a light bending formula 
using a fully relativistic approach, following the route de- 
scribed by Weinberg [T^l . We also develop a formulation 
for extended distributions of sources, which is important 
not only for lensing, but also in future, planned work that 
includes astronomical iV-body simulations using MOG. 

We begin in Section |TT1 reviewing the basic equations 
of MOG and the results of an exact numerical solution 
in the spherically symmetric case. We proceed by devel- 
oping the equations of motion in the vicinity of a MOG 
point source in Section Hill In Section [TVl we obtain an 



exact treatment for light bending in the field of a point 
source. We generalize our discussion to extended sources 
and lensing in Section |Vl Lastly, in Section IVll we con- 
clude by exploring the possible consequences for the Bul- 
let Cluster and for the merging clusters of galaxies Abell 

520 [nl. 



II. MODIFIED GRAVITY THEORY 



Our modified gravity theory is based on postulating 
the existence of a massive vector field, 0^ . The choice of 
a massive vector field is motivated by our desire to in- 
troduce a repulsive modification of the law of gravitation 
at short range. The vector field is coupled universally 
to matter. The theory, therefore, has three constants: 
in addition to the gravitational constant G, we must also 
consider the coupling constant lo that determines the cou- 
pling strength between the (j)^ field and matter, and a 
further constant n that arises as a result of considering a 
vector field of non-zero mass, and controls the coupling 
range. The theory promotes G, /z, and uj to scalar fields, 
hence they are allowed to run, resulting in the following 
action P, Q]: 



where 
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where Sm is the "matter" action, i?^^ = 9^10^ — d^(j)f_i, 
while Vip{(j)), Vc{G), and y^(^) denote the self- 

interaction potentials associated with the vector field and 
the three scalar fields. The symbol is used to denote 
covariant differentiation with respect to the metric g^'', 
while the symbols R, A, and g represent the Ricci-scalar, 
the cosmological constant, and the determinant of the 
metric tensor, respectively. We define the Ricci tensor as 

R^ii^ ~ da^"i, — di.T'^^ + r^^r^^ — r^^r^j,. (5) 

Unless otherwise noted, our units are such that the speed 
of light, c = 1; we use the metric signature (+, — , — , — ). 

In the case of a spherically symmetric field in vacuum 
around a compact (point) source, we were able to derive 
an exact numerical solution [2] . We found that the scalar 
fields G, fi, and uj remain constant except in the immedi- 
ate vicinity of the source. The spatial part of the vector 
field (j)f_i, is zero, while its t-component obeys a simple 
exponential relationship. Meanwhile, the metric is ap- 
proximately the Reissner-Nordstrom metric of a charged 
source. 

Specifically, given a spherically symmetric, static met- 
ric in the standard form 



'■{de"^ +sm^9d(P^), (6) 



we found that, for a source mass M, 
A ~ B-\ 
1 



2GnM AttloGoQI 



Go — Gn + (Goc — Gn)- 
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B ~ -----o ^ (8) 

G ~ Go — Gn + (Goc — Gn)—^= (9) 

(10) 
(11) 

V 

0t ~ -Q5^, (12) 
r 

where Gn is Newton's constant of gravitation, Qs = kM 
is the fifth force charge associated with the source mass 
M, while Goo, D and E are constants. Further, 
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When r is large (that is, large relative to the 
Schwarzschild-radius rg — 2GnM for a source mass M), 
the metric coefficients become 
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B 
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III. EQUATIONS OF MOTION IN A 
SPHERICALLY SYMMETRIC FIELD 



To develop an equation of motion for a point particle, 
and use it to derive a formula for light bending, we follow 
the approach presented in [lOj . 

We begin with the point particle action in MOG, which 
is written in the form 



S = — J {m + au!q5(t>i^iU^) dr 

= - / (m^ gapu°-uP -t- auq^cfif^u^^) dr, (19) 



where m is the point particle mass, and is its fifth force 
charge. The fifth force charge is assumed to be propor- 
tional to m, such that — nm, where k is a constant. 
In earlier work, we determined that k = y^GN/to- To de- 
velop an equation of motion, we compute the derivatives 
of the Lagrangian with respect to positions and veloci- 
ties: 



(20) 



dC 1 



dx 



-— = -mgaiyU - aojqsipiy, 



(21) 
(22) 



d dC du"^ o 

-TTTr-r = -mgau— mu'^gau.pu 

at ou" dr 

-aujq5U^4>^jj - aujjjq5u'^(l)i,. (23) 
We can now construct the Euler-Lagrange equation: 



dC _ d dC 
dx'^ dr du'^ 



mga 



du" 



+ mVapuU'^u'^ + at^q5u'^{4>y,(3 - 4>l3,iy) " CtUl ^^q^cj) ^^u'^ + UUl ^aqb(t>vU°' = 0, (24) 
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or, after rearranging terms, 



where we used Bap = 4>p a — 4'a p- This is the same as 

Eqs. (31-32) in ' ' R' R' 

If uj is constant (as confirmed by the numerical solution -p^ _ r* — -— 

in 0), we get ' "~2A' '-"25' ^''^ 

^+r>"uM=ag5C.«V"i?a,. (26) 

^ no^-j, r^^ = cot0, (29) 

The Christoffel-symbols associated with the spherically sii^2 ^ ^ 

symmetric metric ^ are = -^^ — , F^^ = — cos 6* sin 6'. (30) 

The equations of motion read, using = nm and di- 
viding through by m, 



eft B' dt dr 3 



fdty A' (dry r sin' 6 fd4>Y 

[-d-r) ^ 2A[-d-r) - ^ [-d-r) = --^-'^B.,. (32) 
d-'e ^2drde_^^^^^^Jdf^' ^ (33) 



dr^ r dr dr \ dr 

d^ct) 2drd(t) 2cote'^^ 
dr'' r dr dr dr dr 



d'^d) 2 dr dch dO deb a aj, , , 

^ ■ ^ ' 2cot6l— -- aKUJul^g't"f'B^p. (34) 



We can set 6 ~ 7r/2 without loss of generality. We can also recognize that the only non-zero components of Bap are 
Btr = ~Brt = dt(t)r ^ 9^0* = ~4>'f Wc get 



d^t B' dt dr anuj dr 



dr'^ B drdr B dr 

d'^r B' f dtV A' fdrV r ( dd>\ aKLU dt ,, 



dr"^ 2A \ dr 2A \ dr A \ dr A dr 



(35) 

t, (36) 



^ + -?^-0. (37) 
dr^ r dr dr 



Eq. (1551) '^^^ be rearranged after multiplying with B: 
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dB dt 
dr dr 
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dr'^ dr dr dr 
which can be integrated to yield 

C — aKOjipt 



dt 
d^ 



B 



(38) 



(39) 



Since B ^ 1 and 



as r ^ oo, an asymptotically 



flat spacetime requires C = 1. 

Eq. ((37)) can be integrated directly: 



r^^ = J, 



dr 



(40) 



where J is a constant of integration, which we identify 
as the angular momentum per unit mass. 

Using these results in Eq. (|55)) yields 



In the case of general relativity, the right-hand side of this 
equation is zero, and the particle moves along a geodesic. 
This is not the case here, but we can still integrate our 
equation. Multiplication by 2Adr/dT leads to 



J2 (1 - aKUJ(t)tf 



B 



= 0. (42) 



Integration yields 



drj ■ B - (4^) 

where £ is another constant of integration. After using 



, we get 
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'dtj 



J2 



(1 — aKU!(j)t)^ 



1 

B 



£ 



(1 — anuxptY 
(44) 



This is an exact result. From this result, we can de- 
velop the equation of motion for a non-relativistic particle 
in the usual form. For this, let us assume that we are far 
from a source, and the metric is that of Schwarzschild, 
in accordance with (fT7| and psp . Then, our equation of 
motion becomes 



1 fdrV J2 1 1 -£ 



(1 - 2GM/rf V di) ^ (1 - aKtj0t)2 ^2 1 - 2GM/r (1 - aKuj(t)tY ' 
Multiplying both sides with (1 — 2GM/r)^ gives 

dt) ^ (1 - aKUj<Pt)r'^ \ ~) ~\ ~) " (1 - aKUj(f>tf \ ~) 

Fully differentiating with respect to t, dividing through with 2dr/dt and then rearranging terms yields 

d^r (1 - 2GM/rf (1 - 2GM/rf (/)' + ~ 2GM/rf SJ^GM / 2GM\ 2GM 

-£{\-2GM/rf ,, ?,£{l-2GM/rf GM 

= ~ri Tu""'*'^'^* — n Tv2 —■ (47) 

(1 — anujqjt) — 



From taking the large-r limit, we get 

£ = l-v^, (48) 

where v = dr/dr ~ dr/dt is the velocity of the particle 
at infinity. For photons, £" = 0, for material particles 
in unbound orbits, < £ < 1, and for bound particles, 
£>l. 

For a non-relativistic particle, ^ 1. Further, in the 
weak field limit, 1 — 2GM/r ~ 1, 1 — aKuj(pt ~ 1, and 
j2/2r2 < 1, so pll| becomes 



d^r J2 . GM 



(49) 



IV. THE BENDING OF LIGHT 

To calculate light bending, we return to the exact re- 
sult in ([331). The first step is eliminating dr using (|40p . 
since we are interested in the shape of the orbit, not its 
time evolution: 



A 



dr 



' dr 



J2 (l-aKW0f)2 
72 



B 



or, 



A dr 



£ (1 — aKuxjjtY 

j2 + 7^ • 
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(51) 
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From this we get 



vanishes. Then, ([51]) becomes 



^1/2 



T7J dr. (52) 



[1 - aKUj(t)t{ro)]^ 
J = ro\l hv^ - 1. 



B 



At closest approach to a source, r = tq and dr/dcj) Putting this into ([52|) gives 
I 



(53) 



^1/2 



-,1/2 



dr. 



The deflection for a particle coming from infinity to tq and then off to infinity is twice this angle: 
For a photon, v ~ I and 



Acj)^ = 2 



^1/2 



1 f [l-aKi^0t]^ 1 f [l-Q;Ka;0t(ro)]- 
^ I B / I Biro) 

In the case of weak fields, 1 — aKLU(j)t ~ 1 and we get 



-,1/2 



rl AB A 



^1/2 



dr 



(54) 



(55) 



(56) 



(57) 



This formula is formally identical to the light bending formula in the weak field limit of general relativity, with one 
notable difference: instead of G = Gn, we are using G = Goo = (1 + C()Gn in the Schwarzschild coefhcients A and 
B. From this formula, the approximate deflection can be calculated as P. [Toj: 



A(f> = 



4GM _ 4(1 + a)GNM 



ro 



ro 



If (j)t cannot be ignored, we can use the form 0]: 



-fir 



-fir 



= -nM- 



(58) 



(59) 



where — kAI is the fifth force charge of the source with mass M. This yields the formula for light bending in the 
strong field of a point source in the form 



A(j)^ = 2 
where we used k^lu = Gat. 



[r + aGArMexp(-/ir)]2 B{ro) 1 
[ro + aGwM exp(-/xro)]2 AB A 



-1/2 



dr 



(60) 



V. EXTENDED SOURCE DISTRIBUTIONS 



such that J — 0: 



The nonrelativistic equation of motion (|^^ can be fur- 
ther simplified when only radial motion is considered. 



GnM 



[l + a - q;(1 + /ir)e 



(61) 
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This corresponds to the potential 



-fir 



(62) 



such that (r = |r|): 

r = -V^W. (63) 
We write the potential as the sum of two constituents: 
$(r) = $^(r) + $y(r), (64) 

where 

{l + a)GNM 



-I>Ar(r) = -- 



and 



$y(r) 



uGnM 



-fir 



(65) 



(66) 



We find that <I>Ar is the solution of the Poisson equation: 

V2$Ar(r) = 4^(l + a)G^vp(r), (67) 

with p(r) = MS^{r), where S is Dirac's delta function. 

For $y, we consider the inhomogeneous Helmholtz 
equation: 



This equation is solved by 



G(r) 



47rr 



(68) 



(69) 



Replacing k with z/i, we find that 

(V^ - A/^)$y(r) = -AiraGNPir). (70) 
Adding ^ and we get 

V2$(r) =47rGjvp(r)+A*^$y(r). (71) 

If /i = 0, we get back the Poisson equation for Newtonian 
gravity, as expected. 

If p is not a point source but a general continuous 
distribution of matter, the solution of the inhomogeneous 
Helmholtz equation is in the form 

$y(r) = aGN J y^P(f) d^^- (72) 



Thus, 



V2$(r) = 4^^GNp{r)+a^i^G 



N 



r — r 



— P(f)rf f. (73) 



Or, since r = — V$(r), we can write 



V • r = -AnGNPir) - ap.'^G 



N 



-/j|r-r| 



— p(f) d-^r. (74) 



FIG. 1: The position r of a test particle relative to a point x 
inside an extended mass distribution characterized by p(x). 



In the case of a point source with mass M, fi is given 
by (Uni), while a = {Gq - Gn)/Gn is determined by 



iVM + E)^ \Gn 



1 



(75) 



Clearly, in the case of an extended mass distribution, 
these expressions must be modified. Presently, we do not 
have solutions to the MOG field equations for extended 
matter distributions. Therefore, we treat this problem 
phenomcnologically. We seek an effective mass function 
7W(x, r), to be used in place of M in ([75)) and (fTO)) . that 
determines an "effective mass", to be used in the for- 
mulae for a and /x, when calculating the gravitational 
influence of matter at x on a test particle located at r 
(see Figure [1]). The function should yield the mass of the 
source in the case of a point source, and a mass propor- 
tional to volume in the case of a constant distribution. 
One function that satisfies these criteria is in the form, 



M (x, r) = / p(x) exp ( -^r^ ) ^'^^ (76) 



with the coefficient ^ to be determined from observation, 
e.g., by comparison with the Bullet Cluster data. For 
a constant distribution p(x) = po, this function gives 
A^(x, r) oc |r — xp. If p is not constant, we get the 
following expressions for a and p.: 

I \ Al(x,r) (Goo A 



//(x, r) 



D 



(78) 



This means that in (j73p . we must now move a and p 
under the integral sign: 



V2$(r) = 47rGArp(r) 



(79) 



g-M(f,r)|r-r| 

+ Gat / a(f,r)p2(f,r)— — p{r) d^r. 

r — r 



Given that A4 depends only on x and ro = |x — r|, we 
can also write (1761) in the form 



X(x,ro)= y p(x)exp(^-e^^^j d'*, (80) 
and a and p can also be written as functions of x and tq . 



D, = D, + D„ 



FIG. 2: The geometry of thin lensing: a foreground object, at distance Di from the observer, deflects hght from a background 
object, at distance Ds, by an angle A(/>. 



For lensing, we can use the K-convergence formula [IS 



4^G(x) DiDi, 
c2 



■p(x) dz, 



(81) 



where Di is the distance from the observer to the lensing 
plane, Ds is the distance from the observer to the distant 
light source, Dis = Ds — Di (see Figure^, and a rectihn- 
ear coordinate system is used such that x = (x, y, z), and 
the 2;-axis coincides with the line connecting the observer 
with the distant light source. In our case, G is given by 



G(x) = [l + a(x,ro)]GAr, 



(82) 



where now represents the distance between the point 
where a straight line connecting the light source with the 
observer intersects the lensing plane, and the point where 
the bent ray of light intersects the lensing plane. For a 
point mass, this distance can be calculated as 



ro 



■ sm Acp ~ A0. 



Ds 



Ds 



(83) 



In the case of extended distributions, A(f> can be calcu- 
lated by 



J V_L«'(x) dz 



(84) 



where Vj^ represents the gradient operator in the (x, y) 
lensing plane. *&(x) would, in the case of general rel- 
ativity, be given by the Poisson equation V^<i>(x) = 
4:TtGnp{x). Our result for light bending (f58|) makes it 
clear that in the case of MOG, we cannot use the nonrel- 
ativistic equation (|7ip . Instead, the correct equation to 
use for the effective potential $ is 



V'$(x) = 4^[1 + a(x, ro)]GNp{^), 



, we get 



or, after substituting ([851) 
V2$(x) 47rGjvp(x) 



1 + a X 



DiDis 2 
Ds c2 



J Vi$(x) dz 



(85) 



(86) 



an equation that is solvable for 

VI. CONCLUSIONS 

We derived an exact formulation for light bending in 
the spherically symmetric field of a point source in modi- 
fied gravity (MOG). Introducing a phenomenological ap- 
proach to model extended distributions of matter, we 
showed how the theory can be used to compute lensing. 
This work is directly applicable to the case of the Bullet 
Cluster 1E0657-558. 

The merging of three clusters of galaxies in Abell 520 
has created a mystery regarding the nature of dark mat- 
ter and its role in merging clusters jTlJ. In contrast to 
the Bullet Cluster, the data show that the dark matter 
does not separate from the baryon plasma situated at ap- 
proximately the center of the merging clusters, whereas 
the galaxies and stars have separated to the sides of the 
cluster. Moreover, these galaxies do not reveal the exis- 
tence of dark matter haloes. In our derivation of the n- 
convergence lensing, we believe that MOG can provide a 
natural explanation for these results without exotic dark 
matter, by modifying the mass profile at the center of 
the merging clusters. This will be the subject of future 
research. 

The phenomenological formalism we developed for 
modeling extended distributions of matter may also be 
useful in future work, including A'^-body simulations us- 
ing MOG, and cosmological computations. 
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